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Abstract 

We study the spontaneous excitation of a uniformly accelerated two-level atom non-linearly 
coupled to vacuum Dirac field fluctuations using the formalism proposed by Dalibard, Dupont- 
Roc and Cohen- Tannoudji (DDC) and generalized by us to the present case in the current paper. 
We find that a cross term involving both vacuum fluctuations and radiation reaction appears, 
which is absent in the linear coupling cases such as an atom interacting with vacuum scalar or 
electromagnetic fluctuations. Furthermore, the contribution of this term actually dominates over 
that of radiation reaction. Thus, the mean rate of change of the atomic energy can no longer be 
distinctively separated into only the contributions of vacuum fluctuations and radiation reaction 
as in the scalar and electromagnetic cases where the coupling is linear. Our result shows that a 
uniformly accelerated atom interacting with vacuum Dirac fluctuations would spontaneously excite 
and a unique feature in sharp contrast to the scalar and electromagnetic cases is the appearance 
of a term in the excitation rate which is proportional to the quartic acceleration. 
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I. INTRODUCTION 



Spontaneous excitation is one of the most interesting and prominent phenomena in the 
interactions of atoms with radiation and so far mechanisms such as vacuum fluctuations 
[1, 2], radiation reaction [3] or a combination of them [4, 5] have been proposed as its 
possible physical explanation. The ambiguity arises as a result of the fact that there exists 
a freedom in the choice of the ordering of the atom and field variables in a Heisenberg 
approach to the problem. The controversy was resolved by Dalibard, Dupont-Roc and 
Cohen- Tannoudji (DDC) [6, 7] who showed that when a symmetric operator ordering is 
chosen then the contributions of vacuum fluctuations and radiation reaction to an atomic 
observable can be distinctively separated when a linear coupling between the atom and the 
field is considered, and furthermore this separation makes them separately Hermitian. Using 
the DDC prescription, one can show that for inertial ground-state atoms, the contributions 
of vacuum fluctuations and radiation reaction to the rate of change of the mean excitation 
energy cancel exactly and this cancellation forbids any transitions from the ground state 
and thus ensures atom's stability in vacuum. While for any initial excited state, the rate 
of change of atomic energy acquires equal contributions from vacuum fluctuations and from 
radiation reaction [8]. The DDC formalism was then fruitfully generalized, in the case of 
linear coupling, to study the radiative properties of atoms in non-inertial motion [9-19] or 
in a thermal bath [20, 21] in a flat spacetime and static ones in a gravitational field which 
interact with vacuum quantum scalar and electromagnetic fields [22, 23]. These studies have 
shed some light on our understanding of the Unruh effect and Hawking radiation from a 
different physical perspective. 

In this paper, we plan to study the spontaneous excitation of a uniformly accelerated 
atom interacting nonlinearly with fluctuating vacuum Dirac fields. Our interest in this issue 
is two fold. First, in the case of an atom in interaction with Dirac fields, the simplest 
Lorentz scalar interaction Hamiltonian that we can introduce is non-linear (as that will be 
given later), and we want to see what happens to DDC formalism when linear couplings 
are replaced by a non-linear one, e.g, will the contribution to the rate of change of the 
mean atomic energy still be distinctively separated into vacuum fluctuations and radiation 
reaction? Second, when we go from the scalar field to the electromagnetic field, we get extra 
contribution proportional to the acceleration squared besides the thermal term characterized 
by the Plackian factor [8, 9, 11, 12]. It is then natural to wonder whether new features will 
appear for a spin 1/2 Dirac field? 

The paper is organized as follows. In Sec. II, we generalize the DDC formalism to Dirac 
fields where the coupling between an atom and the field is non-linear and then, in Sec. Ill, 
we use the generalized DDC formalism to calculate the spontaneous excitation rate of a 
uniformly accelerated atom with fluctuating Dirac fields in vacuum. We conclude in Sec. 
IV. Natural units h — c — 1 and with metric signature (+, — , — , — ) will be used throughout 
the paper. 
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II. GENERAL FORMALISM 



Consider a two-level atom in interaction with vacuum Dirac field fluctuations in four 
dimensional Minkowski spacetime. Let |— ), |+) denote the atomic ground state and excited 
one with energies being —\ojq and +\ojq respectively and rc M = (x°,x) = (t,x,y,z) denote 
the Minkowski coordinates referring to an inertial reference frame. Assume that the atom 
is on a stationary trajectory x(t) = (t(r),x(r)) where r represents the proper time. 

In the Dicke's notation [24], the Hamiltonian that governs the evolution of the atom with 
respect to r is 

H A (r) = cu R 3 (t) (1) 

where #3(0) = ||+)(+| — \\— )(— |- The field ip(x) that the atom is assumed to be coupled 
to satisfies the Dirac equation 

(% - m)^(x) = (2) 

where ^ = 7 M <9 M with = and 







' aA 








[0 -I, 




[~<Ji J 



7°= , 7 f = (3) 

\0 -I J \-<7i J 

with a being the Pauli matrices. The 7 matrices satisfy the algebra: {7 M ,7 1/ } = 2g^ u where 
{ , } denotes the anti-commutator. Solving the Dirac equation, we can expand the field 
operator in terms of a complete set of plane-wave solutions as 

d 3 k 



in which b(k, s, t) and d + (k, s, t) are respectively annihilation and creation operators of par- 
ticles and antiparticles with momentum k and spin s and 

u(k,s) = J + m = u(0, s), (5) 
y/2m{uj: + m) 

v(k,s) = ^ + m - v(0,s) (6) 

where jt = fc M 7 M , u(0,s) and v(0,s) are the unit spinors in the particle rest frame with 
k = (m, 0) which have only upper or only lower components. When there is no coupling 
between the atom and the field, b(k, s, t) = b(k, s)e~ wt ) d(k, s, t) = d(k, s)e~ VjJt . The vacuum 
is then defined by the annihilation operators as 

6(jfe,s)|0)=d(jfe,s)|0) = 0. (7) 
The annihilation and creation operators satisfy the anti-commuting relations with 

{b(k,s),b + {k',s')} = {d{k, s),d + {k / , s')} = S 3 {k - k')S ss > (8) 



3 



and others being zero. 

The free Hamiltonian that governs the evolution of the free Dirac field with respect to 
the proper time r is given by 

H f (t) = J2jd 3 k uj: [b + (k, s) b(k, s) + d + (k, s) d(k, s)]^ . (9) 

s 

The interaction Hamiltonian that describes the coupling between the atom and the field is 
assumed to be 

Hl (r) = hR*{t)${x{t))iI){x{t)) (10) 

where i/}(x(t)) = i]r '(x(r))7° , \i is the coupling constant that is assumed to be small and 
R 2 (0) = \i[R-(0) - R+(0)} with R+(0) = |+)(-| and R-(0) = |-)(+| being the atomic 
raising and lowering operators respectively. These operators obey the angular momentum 
algebra, i.e., [R 3 ,R±] = ±R±, [R + ,R-] = 2R 3 where [ , ] denotes the commutator. Notice 
that here the Hamiltonian is quadratic in the field operator, so the coupling is non-linear in 
contrast to the case of scalar and electromagnetic fields [8, 11] where the coupling is linear. 
With this kind of nonlinear interaction, atomic transitions can occur both via absorption 
and emission of Dirac particle-antiparticle pairs and inelastic scattering of a particle or 
antiparticle even at the lowest order of perturbation. This is in contrast to the linear coupling 
case, e.g., a scalar field, where the quantum is singly absorbed or emitted and inelastic 
scattering occurs only at higher orders. Let us note that this interaction Hamiltonian has 
been used for the Unruh particle detector which is aimed to detect vacuum Dirac field 
fluctuations and responses of the detector are found [25-27]. 

The total Hamiltonian of the system (atom+field) is composed of the above three parts 

H(r) = H A (r) + H F (r) + Hj(t) . (11) 

Starting from the above Hamiltonian, we can obtain the Heisenberg equations of motion for 
the dynamical variables of the atom and the field: 

±iu R±(r) + ifi[R 2 (T), R ± {t)\^{x{t))^{x{t)) , (12) 

i^R 2 (r),R 3 (r)]i}(x(r))^(x(r)) , (13) 

dr 

-iu^k, s, t(r)) + wR 2 (t)&(x(t))4>(x(t)), b(k, s, t(r))]— . (14) 

The Heisenberg equation of motion of the annihilation operator of antiparticles is obtained 
by replacing b(k, s,t(r)) with d(k, s,t(r)) in Eq. (14). In the above equations, we have kept 
the commutators unevaluated for latter use. 

The solutions of the above equations of motion can be split into two parts: the free part 
that exists even when there is no coupling between the atom and the field and the source part 



dR±(r) 

dr 
dR 3 (r) 

dr 

db(k,s,t(r)) 
dt 
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that induced by the interaction between them and characterized by the coupling constant 

R ± (t) = R f ± (r) + R s ± (r), 
R 3 (t) = R((t) + R s 3 (r) , 
b(k,s,t(r)) = bf(k,s,t(T)) + b s (k,s,t(r)) . 

Then formal integration of Eqs. (12)-(14) gives 





= R{{T )e ±l ^ T - T °^ , 


(15) 


R s ± (r) 


= in f dr' [R((t'),R{(t)W(x(t'))^(x(t')) , 


(16) 


R((t) 


= Rl(ro) , 


(17) 


R S s(r) 


= in [\t'[rI(t'),rUtW(x(t'))^(x(t')) , 

J T() 


(18) 


b f (k,s,t(r)) 


= b f {k,s,t{T ))e- l ^ T) ' t{T ^ , 


(19) 


b s (k,s,t(r)) 


= in [\r'RUT'W(x(T'))^(x(T')),bf(lsAT))} . 

J TO 


(20) 



In the source parts of the above solutions, all operators on the right-hand side have been 
replaced by their free parts as we have made them only accurate to the first order in the 
coupling constant \i. Similarly, we can derive the corresponding free part and source part of 
the annihilation operators of antiparticles. So, the field operator can be expressed as, 



4>(x(t)) = ^{x{t)) + iJ s (x(r)) (21) 



with 



^(x(r)) = J2j 



d 3 k pfn 
(2^V^ 

x [b f (k, s)u(k, s ) e -<WjE(*W-*(t)))-H*-(*W-*(i>)) 



+rf /t (fc,s)t;(fc,s)e lw s ( ' (r) ^ (To)) - lfe - (i;(r) ^ (To)) ] , (22) 
r(x(r)) = ^ f dr'R{{T')W{x{T')W{x{T'))^{x{T))\ . (23) 



TO 



Our aim now is to study the roles played by vacuum fluctuations associated with ^{x{t)) 
and radiation reaction associated with i/j s (x(t)) in the evolution of the atom. So, we need 
to plug the decomposition (21) into the Heisenberg equation of motion of the atomic Hamil- 
tonian 

= Wi? 2 (r), R 3 (rMx(r))i;(x(T)) . (24) 

In so doing, an issue of operator ordering arises as the free part ^ (x(t)) and source part 
ip s (x{r)) no longer separately commute with the atomic operators. Actually, we can write 

dH A (r) (dH A {r)\ | ( dH A {r) \ | f dH A {r) \ ^ ^ 



dr V dr / , V dr /„_ V dr 

' vj \ / cross 
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where 



(dHAir) 
V dr 

dH A {r) \ 
dr J 



cross 



vf 



itiUJ (\i> f (x(T))4> f (x(T))[R 2 (T),R 3 (T)} 



+ (1 - X)[R 2 (r), R,{r)W {x{r)W {x{r))) , 
if iu (\^(x(T))r(x(r))[R 2 (r),R 3 (r)] 



(26) 



+ (1 - \)[R 2 (t), R 3 (r)W(x(r))r(x(r)) 

+Xr(x(T))^(x(T))[R 2 (T),R 3 (r)} 

+ (1 - \)[R 2 (r), R 3 (T)]r(x(T))4>f(x(T))) 



(27) 



and 



dH A (r) 
dr 



rr 



i^ Q (Xr(x(T))r(x(r))[R 2 (T),R 3 (r)} 



+ (1 - \)[R 2 (r), R 3 (T)}r(x(rW(x(r))) . 



(28) 



Classically, for any A ranging from zero to one, the above equations are equivalent, but quan- 
tum mechanically they are not. We will chose A = 1/2, since such a choice ensures that the 
above three parts are all separately Hermitian. A few comments are now in order. First, the 
rhs of the first and the last equations only involve the free part and source part respectively 
and so they can be interpreted as the sole contribution of vacuum fluctuations and that of 
the radiation reaction respectively. Second, even if we chose a symmetric ordering, there still 
exists a contribution to the rate of change of the atomic energy that involves both the free 
part and source part and this cross term can neither be understood as the sole contribution 
of vacuum fluctuations nor that of the radiation reaction, but rather a combination of them. 
This is in sharp contrast to the cases of linear couplings where no cross term appears and 
the contribution to the rate of change of the atomic energy can be distinctively separated 
into only vacuum fluctuations and radiation reaction [8, 12]. Furthermore, as we will see 
later, the contribution of this cross term dominates over that of radiation reaction. In fact, 
the appearance of a cross-terms is a direct result of the nonlinear coupling between the atom 
and the field. 

Taking the average value of the rate of change of the atomic energy, Eq. (25), over the 
state of the system, |0, b), where represents the vacuum state of the field and b the state 
of the atom, we find that both { dH £^ ) v f and ( dH ^ T ^ ) cross are of the order of fj, 2 , whereas 
( dH ^ T ^ )rr is of /i 3 . So, to the order /i 2 , we have 




dH A (r) 
dr 




dH A (r) 
dr 




(29) 
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with 



dHA(T) \ U^(^{x{t)W{x{t))[R 2 {t),R 3 {t)] + [R 2 (t),R 3 (t)]^(x(t))^(x(t))) , 

(30) 



dr / vf 2 



dr 2 

' cross 



dHA ^\ = ii!,uj (^(x(rW(x(r))[R 2 (r),R 3 (T)} + [R 2 (r), R 3 (r)W (x(r))^(x(r)) 

+r(x(r))^(x(r))[R 2 (r),R 3 (r)\ + [R 2 (r), R 3 (r)}r(x(r))^ (x(r))} . 

(31) 

Here ( ) represents the expectation value over the state of the system. A distinct feature 
as compared to the scalar and electromagnetic cases where the couplings are linear [8, 12] 
is that now the contribution of radiation reaction is of the order higher than that of the 
vacuum fluctuations and thus negligible. Further simplifications by using Eqs. (15)-(18), 
(22) and (23) yield 



d_HAj) 
dr 

dH A (r] 
dr 



= 2^ 2 f T dr'C F (x(T),x(r'))^ X A (r,r') , 



(32) 



2iy? J\t' x F( x ( t ),x(t'))±C a (t,t') (33) 



with 



C f (x(t),x(t')) = ^(Ol^^r))^^^)),^^^))^^))}^} , (34) 

X F (x(t),x(t')) = -^(OK^^Cr))^^^))^^)),^^))] 

+ $ f (x(T'))^ f (x(T')),^(x(T))}^(x(r)))\0) (35) 

being the two statistical functions of the field and 

C A (r,r f ) = l -(b\{R{{r)M{r')}\b) , (36) 
X A (r,r') = l -(b\[R{{r),R{{r')]\b) (37) 
being the two susceptibility functions of the atom which can be further simplified to be 

C a (t, t') = ( fo l^2(0)|rf)| 2 (e lWM(T - T,) + e ~ lu ^ T ' T,) ) , (38) 

d 

X A (r,r') = i^|(6|i? 2 (0)|rf)| 2 (e— -e— ')) . (39) 
d 

The summation in the above two equations extends over the complete set of the states of 
the atom. 
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III. SPONTANEOUS EXCITATION OF A UNIFORMLY ACCELERATED ATOM 



Assume that a two- level atom is initially in state \b) and is uniformly accelerated along 
the stationary trajectory 

t(r) = -sinh(ar) , x(t) = - cosh(ar) , y(r) = z(t) = (40) 
a a 

with acceleration a in four dimensional Minkowski vacuum with quantum Dirac field fluctu- 
ations. An observer moving along such a trajectory is usually called the Rindler observer. 

To evaluate the rate of change of the atomic energy by the formalism generalized in the 
preceding section, we should first calculate two statistical functions of the Dirac field. For 
this purpose, let us introduce the following matrix with respect to the two point function of 
the Dirac field as it would be useful in latter calculations 

S:(x(t),x(t')) = (0\^x(rmx(r f ))\0) 

J ( 27r ) W fc , 

Generically, the above matrix is not stationary even if the world line was. The reason is as 
follows, the spinor, ifj(x(r)), is not "intrinsic" to the observer moving along the world line. It 
will in general "rotate" with respect to the observer's proper reference frame [26, 28], i.e., the 
reference frame Fermi- Walker transported by the observer. To keep the spinor, as it were, 
"in the same direction" with respect to the proper frame reference, it is useful to introduce 
a transformation to the Dirac field operator [28, 29] 

iP(x(t)) -> SW(x(r)) , (42) 
4,{x{t)) ${x{t))S? , (43) 

where the matrix S T is to take care of the Fermi- Walker transportation and it is given, in 
the case of a uniformly accelerated worldline with proper acceleration a, by 

S T = e h a ^i = cos h(ar/2) + 7o 7i sinh(ar/2) . (44) 

Thus it is more appropriate to define a new matrix [26, 28] 

g(T,T') = S T S+(x(T),x(T'))S;, 1 (45) 

which will be later shown to be a function of the time interval Ar = t — t' and can effectively 
simplify the calculations concerned. Similarly, we can define another matrix with respect to 
the two point function of the field as follows 

(S-(x(r'),x(r))) ab = (0\Mx(r))Mx(r'))\0) 

= /^^£* 6 (M)t^ (46) 
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Noticing the relations that 

J2u(k,s)u(k,s) = £±p, (47) 

s 

J2v(k,s)v(k,s) = (48) 



we find 



where 



S+(x(t),x(t')) = (i@ + m)G + (x(T),x(T')) , (49) 
S-(x(t),x(t')) = -(ifi + m)G-(x(T),x(T')) , (50) 

G + {x{r),x{r')) = / — e -<« S (*(r)-t(r'))-H^(*(r)-*(r')) (51) 

and G~(:e(t), rr(r')) = G + (x(t'), rr(r)) are just Wightman functions of the scalar field in four 
dimensional Minkowski spacetime. We can deduce from these relations that them atrices 
S+ and S~ are related as 

S~(x(t),x(t')) = S+(x(T'),x(T))\ m -+- m . (52) 

By using the properties S*" 1 = S- T , S T S T r = S T+T r, and joS T = 5 , _ r 7 , it can be proved that 
the function g(r, r') can be expressed in terms of the interval of the atomic proper time, 
At = t — t', as 

<?(t,t') = g(Ar) 

= (- 1 °d z + mS AT )G + (z(Ar)) (53) 

with 



z(At) = i-sinh^-Ar - ie J . (55) 
For the case of massless Dirac fields, 

G + (z(Ar)) = - 9r t - . (56) 
v v 77 47r 2 [2(Ar)] 2 v 1 

We can show, after lengthy simplifications, that the following functions related with the 
two statistical functions of the Dirac field can be expressed in terms of the above two matrices 

as 

(0|^(x(r))V(x(r))^(x(r / ))V(^(r , ))|0) = Tr[S+(x(r), x(t'))S;(x(t'), x(r))\ , (57) 
(0|[^(x(r'))^(x(O)^(x(r))]^(x(r))|0> = Tr[5 n (x(r),,(r'))S;(x(rV(r))] 

+ Tr[5-(x(r),o;(r'))5-(x(r'),x(r))], (58) 
(0|^(x(r))[^(x(/))^(x(r'))^(x(r))]|0> = -Tr[S+(x(r), x(t'))S;(x(t') : x(r))} 

-Tr[S-(x(r'),x(r))S-(x(r),x(r'))} (59) 
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where TV[---] represents the trace of a matrix. In obtaining the above three equations, 
we have used (7 S T ) 2 = 1 and [A, BC] = {A,B}C - B{A,C}. In Eq. (57), we have 
discarded an infinite constant as it would have no contributions to the integrals in the 
following calculations. Considering massless Dirac field fluctuations (m = 0), we find 



Tr[S:(x(r),x(r'))S-(x(r'),x(r))] = 4[^G + (z(Ar))] 2 . 



(60) 



Similar simplifications can also be obtained for Eqs. (58) and (59). Consequently, the two 
statistical functions of the Dirac field can be shown to be given by 



C F (x(r),x(r')) 
X F (x(t),x(t>)) 



a 



1287T 4 
,6 



+ 



sinh b (f At - ie) sinh b (f At + ie) 



a 



1287T 4 



1 



1 



sinh b (f At - ie) sinh b (f At + ie) 



(61) 
(62) 



Inserting Eqs. (61) and (39) into Eq. (32) and taking the proper time to be infinite long, 
we can express the contributions of vacuum fluctuations to the mean rate of change of the 
atomic energy as 



dH A ( T ) 

dT 



vf 



//V 

1287T 4 



x 



J2\(b\MO)\d)\ 2 oo bd 



f 

Jo 



dAT 



sinh b (|Ar 



+ 



1 



ie) sinh b (f At + ie) 



(63) 



By exploiting the techniques of contour integration and residue theory, the above integral 
can be calculated out to be 



dH A (r) 

dT 



// 



vf 



1207T 3 



2 , , 6 



U) b >UJ d 



1 + 5^ + 4- 



bd 



bd 



i + 
i + 



e 2nui bd /a _ I 

2 

e 2n\u> bd \/a _ 



(64) 



Obviously, just as that of a uniformly accelerated atom interacting with vacuum scalar 
or electromagnetic fluctuations (see Eq. (56) in Ref. [8] and Eq. (28) in Ref. [12]), the 
contributions of vacuum Dirac fluctuations would raise the atomic energy if the atom is 
initially in its ground state (concerned with the term (ub < Ud)) and diminish its energy if 
the initial state of the atom is the excited state (concerned with the term (uj, > to a))- A 
sharp feature in contrast to both the cases of an atom in interaction with massless scalar 
field and electromagnetic field fluctuations is the appearance of a term proportional to a 4 
which is unique to the Dirac fluctuations. 
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Similarly, we can find, by taking Eqs. (62) and (38) into Eq. (33), the contributions of 
the cross term to the mean rate of change of the atomic energy as 



dH A (r) 
dr 



/i 2 a 6 

1287T 4 



£| (61^(0) |d>| 2 w M 

d 

1 



x / dAr 
'o 



sinh b (f At - ie) sinh b (§ At + ie) 



(e 



(65) 



Further calculations give 
(dH A {r) 



dT 



1207T 3 



E \{b\R 2 (0)\d)M d (l + 5-^ + ^) 



1207T 3 



E IflWOMI 



2 , , 6 



1 + 5^ + 4- 



6(1 



(66) 



It is interesting to note that the contribution of the cross term here always diminishes the 
atomic energy no matter if the initial state of the atom is the ground state or the excited 
state, so the cross term just plays the role as that played by the radiation reaction in the 
scalar and electromagnetic field cases [8, 12]. The presence of it ensures that the inertial 
ground state atoms will be stable in the Minkowski vacuum since now the contribution of 
radiation reaction is of order higher than that of the vacuum fluctuations and thus negligible. 

Adding up Eqs. (64) and (66), we obtain the total mean rate of change of the atomic 
energy 



^dH A (r) 



dT 



tot 



60tt 3 



E KW°)I<*>I 



2, ,6 



1,(1 



LO b >UJ d 



l + 5^r + 4- 



bd 



bd 



1 + 



1 



+ 



60tt 3 



E |(fe|^ 2 (0)|rf)| 2 a; 6 6 d (l + 5^ + 4 



e 2-Ku bd /a _ I 
1 



U bd I e 27V ^bd\/a _ i 



(67) 



The above result reveals that uniformly accelerated atoms in ground state in interaction with 
fluctuating Dirac fields in vacuum would spontaneously excite. This is consistent with the 
response of Dirac particle detector found in Ref. [27]. However, our result can be considered 
as providing a transparent underlying physical mechanism on why such a particle detector 
clicks. A distinct feature is the appearance of the a 4 term in addition to the Planckian 
thermal factor which can be viewed as a result of an ambient thermal bath at the Unruh 
temperature at T = a/2ir. This term is absent in both the scalar and electromagnetic cases. 
For a typical transition frequency of a hydrogen atom, u> ~ 10 16 s -1 , the corrections due to 
the acceleration are negligible when a u ~ 10 24 m/s 2 and it becomes appreciable when 
the acceleration approaches ~ 10 24 m/s 2 . In fact, the contribution of the a 4 term which is 
unique to the current case becomes dominant when a ^> u>. Here, it is also interesting to 
note that the response of an accelerated detector to the Dirac vacuum fluctuations near a 
fluctuating event horizon has recently been studied [30]. 
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The above result reduces in the limit a — > to 



{^) tt =-^^ mR ^ w ^- (68) 

This means that no spontaneous excitation would occur for inertial atoms in their ground 
states since the term associated with Ub < vanishes, as a result of the complete cancel- 
lation of the contribution of vacuum fluctuations and that of the cross term. Without the 
cross term, inertial ground state atoms in vacuum would not be stable. This is in sharp 
contrast to the scalar and electromagnetic cases that it is the cancellation of contribution of 
radiation reaction which is negligible in the present case and that of the vacuum fluctuations 
that ensures the stability of the inertial ground state atoms in vacuum [8, 12]. 



IV. SUMMARY 



We have generalized the DDC formalism to the case of an atom interacting with vacuum 
Dirac field fluctuations where the coupling between the atom and field is non-linear, and 
find that a cross term appears, which involves both vacuum fluctuations and radiation 
reaction and which is absent in the linear coupling cases such as an atom interacting with 
vacuum scalar or electromagnetic fluctuations. Furthermore, the contribution of this term 
actually dominates over that of radiation reaction. Thus, the mean rate of change of the 
atomic energy can no longer be distinctively separated into only the contributions of vacuum 
fluctuations and radiation reaction as in the scalar and electromagnetic cases where the 
coupling is linear. But rather, the evolution of the atom is governed in the leading order 
by vacuum fluctuations and the cross term that is a combined effect of vacuum fluctuations 
and radiation reaction. 

We then calculated the spontaneous excitation of a uniformly accelerated atom in in- 
teraction with fluctuating vacuum Dirac fields. Our result shows that such a uniformly 
accelerated atom would spontaneously excite in vacuum and the excitation rate contains, 
besides what can be viewed as a result of an ambient thermal bath at the Unruh temperature 
at T = a/27r, a term proportional to quartic acceleration which is absent in both the scalar 
and electromagnetic cases. 
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